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ABSTRACT 

We present photometric and spectroscopic data analysis and orbital period study of an early- 
type interacting binary system V382 Cyg by using all the available data. We made a simultaneous 
light and radial velocity curve solution. The derived physical parameters of the primary and 
secondary stellar components are Mi = 27.9(5) M Q , M 2 = 20.8(4) M Q , R x = 9.7(2) R Q , R 2 
= 8.5(2) Rq, log(ii/i©) = 5.152(20) and log(L 2 /£©) = 4.954(19) while the separation of the 
components is a = 23.4 Rq. Newly obtained parameters yield the distance of the system as 
1466(76) pc. Analyses of the mid-eclipse times indicate a period increase of ^ = 4.2(1) x 10~ 7 
days/yr that can be interpreted in terms of the high mass transfer (^ = 6.1(5) x 10~ 6 M /yr) 
from the less massive component to the more massive component. Finally we modelled the 
evolution of the components using nonconservative codes and discussed the obtained results. 
The age of the binary system is estimated as 3.85 Myr. 

Subject headings: Stars: binaries -stars: binaries: close - stars: individual: V382 Cyg - stars: fundamental 
parameters -stars: early-type - stars: evolution 



1. Introduction 

Massive interacting binary systems are the pro- 
genitor candidates of some high energy phenom- 
ena in astrophysics (e.g., supernovae). Solar- type 
close binary systems with convective envelopes 
can evolve in different ways than close massive 
stars. Yakut & Eggleton (2005) modelled late- 
type close binary systems assuming conservative 
and non-conservative cases. The authors proposed 
that a detached configuration can evolve into a 
semi-detached configuration followed by a contact 
configuration. Except late-type systems, we have 
less information on the evolution of the interac- 
tive massive binary systems, since there are a few 
binary systems with well-determined parameters. 
Therefore, systems like V382 Cyg, RY Set (Djura- 
sevic et al. 2008) and V729 Cyg (Yasarsoy et al. 
2013, submitted), etc. are crucial in observational 
studies of close massive binaries. 

Photometric studies of V382 Cyg have been 
published by different researchers. The binary 
system V382 Cyg was discovered as a variable 
star by Morgentoth (1935) and Petrov (1946) and 



they classified the system as an eclipsing binary. 
The photometric UBV light variations of the sys- 
tem was obtained by Landolt (1964) by using 16- 
inch KPNO telescope. Modelling the light and ra- 
dial velocity curves of the system, Bloomer et al. 
(1979) derived an orbital inclination of 87° and Mi 
= 26.7 (5)M Q and M 2 = 18.9 (4)M S . Later pho- 
tometric studies were performed by Cester et al. 
(1978), Mayer et al. (1986), Harries et al. (1997), 
Degirmenci et al. (1999) and Qian et al. (2007). 

Spectroscopic studies of the system were pro- 
vided by Pearce (1952), Popper (1978), Koch et 
al. (1979), Harries et al. (1997), Burkholder et 
al. (1997) and Mayer et al. (2002). Popper 
(1978) studying the He I and He II lines derived 
the semi-amplitude of the radial velocity curve and 
the mass function for the components as Ki = 255 
kms" 1 , K 2 = 360 kms -1 , Mi sin 3 i = 26.7 M 
and M2sin 3 j = 18.9 Mq. Later on Popper and 
Hill (1991) recalculated the parameters as K\ = 
276 kms" 1 and K 2 = 400 kms 
as Mi = 32.6 (1.8)M Q and M 2 = 
using photometric and spectroscopic data Harries 



1 and the masses 
22.9(1.3)M . By 
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et al. (1997) derived the mass of the components 
as Mi = 26.0 (7)M and M 2 = 19.3 (4)M S . Re- 
cently, Mayer et al. (2002) derived the masses 
for the stars as M x = 29.2 M Q and M 2 = 21.2 
Mq. Spectral types of the components were de- 
fined by Pearce (1952) as 06.5+07.5, by Popper 
(1980) as 07.3+07.7, by Harries et al. (1997) as 
07.3+07.7 and by Burkholder et al. (1997) as 
06.5V(f) + 06V(f). They gave distance module 
for the system as 11.5(5) mag. Garmany & Sten- 
cel (1992) classified V382 Cyg as a member of Cyg 
OBI association. 

2. New Observations 

The observations were carried out on 25 nights 
in 2007-2011 with the 40-cm telescope that is 
equipped with the 2048 x 2048 Alta CCD cam- 
era at Ege University Observatory (EUO). The 
exposure times are 15 s for B, 5 s for V, 4 s for 
R and 4 s for I filter. The reduction and analysis 
of the frames have been managed using the IRAF 
packages to subtract the bias and dark then di- 
vide flat field, followed by aperture photometry 
(DIGIPHOT/APPHOT). Comparison and check 
stars are selected as HD 193204 (V=8 m 32, B- 
V=0 m 46), CSC 2684 1450 (V=9 m 60, B-V=l m 12), 
HD 228802 (V=9 m 63, B-V=0 m 44) and HD 193344 
(V=7 nl 60, B-V=-0™06). The magnitudes and 
colors are obtained from TYCHO-2 catalogue 
(H0g et al. 2000). We have studied all the nights 
and each frames separately during the data reduc- 
tion and standard deviations are estimated for B, 
V, R and I bands as 0.018, 0.013, 0.016 and 0.015, 
respectively. 3095, 2861, 2871 and 2772 points are 
obtained in B, V, R and / bands, respectively. 
We list all the observed data in Table [TJ The light 
curves of the system obtained in this study are 
shown in Fig. QJ,. 

We obtained three new minima times through- 
out these new observations. They are collated 
with those published and listed in Table [2] with 
their errors. Using Table [5] we derived a new lin- 
ear ephemeris Eq. ([IJ that we believe to be useful 
for future studies on the system: 

HJD Mini = 24 54344.4972(8) + 1.8855270(2) x E. 

(1) 



Table 1: Observational points of V382 Cyg. Helio- 
centric Julian Date, phase, (Am) and their corre- 
sponding filters are listed. (Am) is the magnitude 
difference between variable (V382 Cyg) and com- 
parison star (HD 193204). The phases were cal- 
culated using the Eq. (1). Table Q] is given in its 
entirety in the electronic edition of this paper. A 
portion of it is shown here for guidance regarding 
its form and content. 



HJD 


Phase 


Am 


Filter 


2454292.33107 


0.3334 


0.331 


B 


2454292.33444 


0.3352 


0.330 


B 


2454292.33540 


0.3357 


0.324 


B 


2454292.33636 


0.3362 


0.326 


B 


2454292.33731 


0.3367 


0.325 


B 


2454292.33826 


0.3372 


0.328 


B 


2454292.33922 


0.3377 


0.315 


B 


2454292.34018 


0.3382 


0.326 


B 


2454292.34113 


0.3387 


0.332 


B 


2454292.34209 


0.3392 


0.337 


B 


2454292.34304 


0.3397 


0.330 


B 


2454292.34400 


0.3403 


0.331 


B 


2454292.34495 


0.3408 


0.338 


B 


2454292.34592 


0.3413 


0.329 


B 


2454292.34688 


0.3418 


0.327 


B 
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3. Eclipse timings and period study 
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Fig. 1. — (a) The observed and the computed light 
and (b) radial velocity curves of V382 Cyg. The 
light curves in B, V, and R bands are altered by 
values of 0.75, 0.50, and 0.25, respectively, for a 
good visibility. See text for details. 



The difference between the observed (O) and 
calculated (C) times of minima in an eclipsing bi- 
nary system can provide us information about any 
orbital period variation. The orbital period study 
of V382 Cyg was provided by Mayer et al. (1986), 
Mayer et al. (1998), Degirmenci et al. (1999) 
and Qian et al. (2007). Using times of minima 
Mayer et al. (1986) deduced the linear elements 
of the system. Degirmenci et al. (1999) taking 
into account only the mass transfer between the 
components derived the mass transfer rate from 
the less massive component to the more massive 
one as 5.0 x 10~ 6 M yr -1 . Qian et al. (2007) 
analyzed the O-C diagram assuming both mass 
transfer (M = 4.3 x 10~ 7 M yr _1 ) and a third 
body. 

Unlike previous studies of V382 Cyg, in this 
study we used only the most accurate (photoelec- 
tric and CCD data) minima times. We solve the 
O-C curve to find a parabolic variation which can 
be ascribed as mass transfer from the less massive 
companion to the more massive one. A total of 65 
data points obtained with photometric/CCD ob- 
servations are used to study the period variation 
of the system. Table [2] shows the data used during 
the analysis. The weighted least-squares method 
is used in order to determine the parameters of the 
upward parabolic variation. The residuals (O-C)j 
indicate a quadratic solution (Fig. [5^). In order to 
estimate the light elements, the differential correc- 
tion method is used. The residuals (O-C)tt show 
a sinusoidal variation (Fig. [5}d). Therefore, we as- 
sumed a sinusoidal variation superposed on an up- 
ward parabola during O-C analysis. A sine-like 
variation in the O-C curve, where both the pri- 
mary and the secondary minima follow the same 
trend, suggests the light time effect via the pres- 
ence of a tertiary component. We used Eq. @ 
during the O-C analysis. 

Mml = T + P E+\%E 2 



(2) 

where T Q is the starting epoch for the primary min- 
imum, E is the integer eclipse cycle number, P a is 
the orbital period of the eclipsing binary and 012, 
i' , e' and u>' are the semi-major axis, inclination, 
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eccentricity and the longitude of the periastron of 
eclipsing pair about the third body and v' denotes 
the true anomaly of the position of the center of 
mass. Time of periastron passage T" and orbital 
period P' are the unknown parameters in Eq.Q 
(see Kalomeni et al. 2007 for details). 

4. Simultaneous solutions of light and ra- 
dial velocity curves 

The BVRI light curves obtained in this study 
are solved simultaneously with the radial veloc- 
ity curve of Harries et al. (1997) and Mayer et 
al. (2002) using Phoebe (Prsa & Zwitter 2005) 
which is based on the WD code (Wilson & Devin- 
ney 1971; Wilson 1979; Wilson 1990). During the 
solution we used weighted light curves that are 
constructed according to the standard deviation 
of each filter. The effective temperature of the 
hot star is chosen according to spectroscopic stud- 
ies. The albedos values A\ and A2, for the hot 
and cooler components are adopted from Rucin- 
ski (1969), whilst the values of the gravity dark- 
ening coefficients g\ and gi are taken from von 
Zeipel (1924). The logarithmic limb-darkening 
law is used with the coefficients adopted from van 
Hamme (1993) for a solar composition star. The 
adjustable photometric parameters are orbital in- 
clination, i, surface potential, f2i,2, temperature of 
the secondary component T2, luminosity L\ and 
the mass ratio q. The center of mass velocity Vo 
and semi-major axis a are also set as free parame- 
ters as well as the time of minimum light, To and 
the orbital period, P. Table [4] summarizes the re- 
sult of the analysis. B,V, R and I light curves and 
the velocity curve computed using the determined 
parameters are shown by solid lines in Fig. [T^. 

During solution we assumed both contact (C) 
and semi-detached (SD) configurations. We em- 
phasized that both solutions have similar results 
(see Table 2}. In Fig. Q] the light curve fits for 
a contact and a semi-detached configurations are 
shown with a dotted and a solid line, respectively. 
For the contact model the filling factor, which is 
expressed by (f2 m — fi)/(f2j_ — ^ ou t) that varies 
from zero to unity from the inner (f2j n ) to the 
outer critical surface (£\mt)' * s calculated to be 
/ = 0.09. New solutions indicate that the system 
has a weak contact degree or that it is a semi- 
detached binary with a near-contact configuration. 



Table 2: Times of minimum light for V382 Cyg. 



HJD* Min 


Rcf 


HJD Min 


Rcf 


HJD Min 


Rcf 


2436814.7725 


1 


2448167.4765 


7 


2449913.4766 


17 


2438987.8298 


1 


2448186.3306 


7 


2449913.4785 


10 


2440385.9310 


1 


2448186.3316 


7 


2449915.3654 


17 


Z4t4UOOD.o 1 Dl 


^ 


9/1 /IS/1/17 A 7/1 T 





9/1/1QQ9Q r;(l7/1 


i 7 


2440387.8173 


1 


2448447.4764 


8 


2449930.4476 


17 


2441129.7698 


1 


2448498.3835 


8 


2449931.3909 


17 


2442651.3844 


2 


2448498.3844 


8 


2449945.5342 


17 


2442659.8678 


3 


2448564.3740 


9 


2449946.4757 


17 


2442940.8071 


4 


2448839.6629 


10 


2449947.4157 


18 


2442961.5424 


3 


2448843.4381 


11 


2449947.4171 


17 


2443366.9333 


5 


2448843.4395 


12 


2450671.4673 


19 


2444445.4450 


6 


2448843.4397 


12 


2451413.4316 


20 


2444757.4930 


6 


2448859.4664 


11 


2451429.4568 


20 


2444842.3538 


6 


2448878.3252 


10 


2452817.2130 


21 


2445598.4426 


5 


2449221.4844 


13 


2453837.2977 


22 


2446274.3971 


5 


2449221.4857 


13 


2455430.5686 


23 


2446325.3087 


5 


2449518.4527 


14 


2455483.361 


24 


2446668.4782 


5 


2449534.4780 


15 


2455811.4399 


25 


2447030.4922 


3 


2449550.5029 


14 


2454344.4913(4) 


20 


2447099.3141 


3 


2449567.4771 


16 


2455003.4938(5) 


20 


2447444.3744 


3 


2449568.4222 


16 


2455897.2295(3) 


26 


2448167.4747 


7 


2449880.4846 


17 







References for Table [2] 1- Landolt(1975) 2- Mayer(1980) 3- 
Mayer et al.(1991) 4- Bloomer, Burke, & Millis (1979) 5- Mayer 
et al.(1986) 6- Andrakakou ct al.(1980) 7- Hubschcr, Agerer, 
& Wunder(1991) 8- Hiibschcr, Agerer, & Wunder(1992) 9- 
Blattler(1992) 10- Mayor et al.(1998) 11- Pctor(1992) 12- 
Hiibschcr, Agerer, & Wundcr(1993) 13- Hubschcr, Agerer, & 
Wunder(1994) 14- Peter (1994) 15- Agerer & Hubscher(1995) 
16- Blaettlcr(1994) 17- Dcgirmenci ct al.(1999) 18- Agerer 
& Hiibscher(1996) 19- Agerer & Hiibscher(1998) 20- Agerer 
& Hiibscher(2001) 21- Nagai(2004) 22- Qian et al.(2007) 23- 
Hiibscher(2011) 24- Paschkc(2010) 25- Zasche et al.(2011) 26- 
prcscnt study. 



Table 4: Simultaneous analyses results of the light 
and radial velocity curve and their formal la er- 
rors for V382 Cyg. The column headers C and 
SD refer to the contact model and semi-detached 



model, respectively. See text for details. 

Parameter C SD 



in 


84.5(1) 


85.3(2) 


q = M h /M a 


0.7439(14) 


0.745(2) 


a (R ) 


23.45(13) 


23.47(12) 


Vo (kms- 1 ) 


7.1(1.6) 


6.3(1.6) 




3.275(3) 


3.390(9) 


n 2 


3.275(3) 


3.320 


Ti (K) 


36000 


36000 


T 2 (K) 


34415(270) 


34578(240) 


Fractional radius of primary comp. 


0.4139(5) 


0.4087(2) 


Fractional radius of secondary comp. 


0.3604(6) 


0.3526(2) 


Ai = A 2 


1.0 


1.0 


91 = 92 


1.0 


1.0 


Luminosity ratio: L (%) 






B 


57.9 


56.2 


V 


58.0 


56.3 


R 


57.7 


56.1 


I 


57.5 


55.8 
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Table 3: Orbital elements of the tertiary component and mass transfer ratio for V382 Cyg. The standard 
errors 1<7, in the last digit are given in parentheses. 



P ar clihg t e r 


Unit 


Value 




[HJD] 


2436814.7599(59) 


Po 


[day] 


1.885515(2) 


P' 


[year] 


43.9(1.7) 


r 


[HJD] 


2419843(1611) 


e! 




0.41(8) 


u' 


[°] 


34(15) 


a\2 sini' 


[AU] 


2.49(14) 


/(rn) 


[M ] 


0.0080(2) 


TO3 ; i' = 30<> 


[M ] 


5.7 


"J3;i'=90° 


[M ] 


2.8 


1 dP 

2 dE 


[c/d] 


1.08(1) x 10~ 9 



5. Results and Conclusion 

Newly obtained B, V, R and I light curves of 
the early-type interacting system V382 Cyg have 
been solved with earlier published spectroscopic 
studies. We derived the orbital and the physical 
parameters of the components from simultaneous 
solution and listed the physical parameters in Ta- 
ble [5] We calculated the masses as 27.9 M Q and 
20.8 M Q which are ~ 7% higher than those given 
by Harries et al. (1997) and ~ 7% lower than 
Mayer et al. (2002) results. The distance of the 
system to the Sun is estimated as 1455 =F 76 pc us- 
ing the observed parameters. During the calcula- 
tions bolometric corrections are taken from Lanz 
& Hubeny (2003), the effective temperature and 
absolute magnitude of the Sun are taken as 5777 K 
and 4.732 mag. Garmany & Stencel (1992) pro- 
posed that V382 Cyg may be a member of the 
Cyg OBI association. The distance to Cyg OBI 
reported by Uyamker et al. (2001) in the range 
from 1.25 kpc to 1.83 kpc. In this study, the dis- 
tance of V382 Cyg was found to be within this 
range, so it is likely to be a member of Cyg OBI. 

We have collected 65 times of mid-eclipses in- 
cluding our three new ones to search any peri- 
odic variation. Analysis of O-C residuals can 
give us the reason of a periodic and/or non- 
periodic variation. We applied the least-square 
method to the Eq. ([2} to find the period change 
rate. Then the mass transfer rate is calculated 
and the possible third body parameters are esti- 
mated (Table [3]). According our analysis, there 
is a mass transfer from the less massive (ini- 



tially massive one) to the more massive component 
= 6.1 x 10~ 6 M Q /yr) and the period change 
rate (£ = 4.5 x 10 6 ). In this study, the obtained 
mass transfer ratio is nearly ten times higher than 
that of Qian et al. (2007). This difference is proba- 
bly due to their formalism. Additionally, we found 
new parameters for the tertiary component with 
43.9 yrs outer period. 

Isolated massive stars and massive stars in bi- 
nary systems (e.g. V382 Cyg) evolve in different 
ways. The evolutionary path towards the final su- 
pernovae will become significantly different. The 
nature of the final stellar remnant may also be 
altered so that while a neutron star might be ex- 
pected, mass transfer may lead to a black hole for- 
mation. Observations of systems similar to V382 
Cyg can allow us to reveal the nature of the bi- 
nary evolution and may help us to understand the 
possible variation of stellar lifecycles (see Eggleton 
2010, Eldridge & Stanway 2009 and Yakut et al. 
2013). Using our newly obtained physical parame- 
ters we provide an evolution model of the massive 
interacting binary system V382 Cyg. We used the 
TWIN version of the EV code (Yakut & Eggleton 
2005, Eggleton 2008, Eggleton 2010) that has been 
developed by Peter P. Eggleton. 

We run a few models using different initial pa- 
rameters. Among these models the best agreement 
with the observations is obtained for the model bi- 
nary system whose initial period is 1.72 days, with 
a primary and secondary masses 28 M a and 23.5 
Mq. The first 28 models shrink slightly since the 
ZAMS is somewhat artificial. The massive compo- 
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Table 5: Astrophysical parameters of the system. The standard errors la in the last digit are given in 
parentheses. 



Parameter 


Unit 


Primary 






Secondary 


Mass (M) 


M 


27.9(5) 






20.8(4) 


Radius (R) 


R Q 


9.7(2) 






8.5(2) 


Temperature (T e g) 


K 


36000 






34415(270) 


Luminosity (L) 


log(L*/L©) 


5.152(20) 






4.954(19) 


Surface gravity (logg) 


cms -2 


3.91 






3.90 


Bolometric magnitude (M&) 


mag 


-8.15 






-7.65 


Absolute magnitude (My) 


mag 


-4.75 






-4.38 


Period change rate ( P) 


d/yr 




4.2(3) x 


io- 7 




Mass transfer ratio (M) 


M Q /yr 




6.1(4) x 


10~ 6 




Distance (d) 


pc 




1466(76) 





nent fills its Roche lobe and RLOF begins at model 
123. At model 201 the parameters are very con- 
sistent with the current observed parameters (Ta- 
ble [6]) . After two more models it reaches a contact 
phase. This configuration is consistent with the O- 
C analysis (see §3). Therefore, our models show 
that the system looks like a semi-detached system 
but very close to a contact phase. We plotted our 
result in the H-R diagram and Mass- Radius planes 
(Fig. [3]). We found the age of the system to be 3.85 
Myr which is in the range of age of Cygnus OBI 
(e.g., Massey et al. 1995). 

Acknowledgments 

The authors would like to thank Peter P. Eggle- 
ton for his valuable comments and suggestions to 
improve the quality of the paper and for providing 
the current version of the EV(TWIN) code and, 
B. Kalomeni for helpful comments. We are very 
grateful to an anonymous referee for his/her com- 
ments and helpful constructive suggestions which 
helped us to improve the paper. This study was 
supported by the Turkish Scientific and Research 
Council (TUBITAK 111T270) and the Ege Uni- 
versity Research Fund. KY acknowledges support 
by the The Turkish Academy of Sciences (TUBA). 
The current study is a part of PhD thesis of B. 
Yasarsoy. 

REFERENCES 

Agerer, F. & Hiibscher, J. 1995, IBVS, 4222, 1 
Agerer, F. & Hiibscher, J. 1996, IBVS, 4383, 1 



Agerer, F. & Huebscher, J. 1998, IBVS, 4606, 1 

Agerer, F. & Hiibscher, J. 2001, IBVS, 5016, 1 

Andrakakou, M., et al. 1980, BBSAG, 49, 1 

Blaettler, V. 1994, BBSAG, 107, 1 

Blattler, E. 1992, BBSAG, 99, 1 

Bloomer, R. H., Burke, E. W. & Millis, R. L. 1979, 
BAAS, 11, 439 

Burkholder, V., Massey, P. & Morrell, N. 1997, 
ApJ, 490, 328 

Cester, B., Fedel, B., Giuricin, G., Mardirossian, 
F. & Mezzctti, M. 1978, A&AS, 33, 91 

Dcgirmenci, O. L., Sezer, C, Demircan, O., Er- 
dem, A., Ozdemir, S., Ak, H. & Albayrak, B. 
1999, A&AS, 134, 327 

Eggleton, P. P. 2008, IAUS, 246, 228 

Eggleton, P. P. 2010, NcwAR, 54, 45 

Eldridge, J. J. 2009, MNRAS, 400, L20 

Garmany, C. D. & Stencel, R. E. 1992, A&AS, 94, 
211 

Harries, T. J., Hilditch, R. W. & Hill, G. 1997, 
MNRAS, 285, 277 

H0g, E., Fabricius, C, Makarov, V. V., et al. 2000, 
A&A, 355, L27 

Hiibscher, J., Agerer, F., & Wundcr, E. 1991, 
BAVSM, 59, 1 



6 



Table 6: Observed and model parameters of the components. 





Star 


M (M ) 


log R 


logT 


log L 


P (days) 


Age (yr) 


Observed 


A 


20.8 


0.9269 


4.5367 


4.954 


1.8855 






B 


27.9 


0.9868 


4.5563 


5.152 






Model No 28 


A 


28.0 


0.8687 


4.5783 


5.0039 


1.6960 


8.6xl0 3 




B 


23.5 


0.8252 


4.5532 


4.8163 






Model No 201 


A 


20.83 


0.9224 


4.5401 


4.9581 


1.8860 


3.85xlO B 




B 


27.65 


0.9741 


4.5606 


5.1440 







Hiibscher, J., Agerer, F., & Wunder, E. 1992, 
BAVSM, 60, 1 

Hiibscher, J., Agerer, F., & Wunder, E. 1993, 
BAVSM, 62, 1 

Hiibscher, J., Agerer, F., & Wunder, E. 1994, 
BAVSM, 68, 1 

Hiibscher, J. 2011, IBVS, 5984, 1 

Kalomeni, B., Yakut, K., Keskin, V., Degirmenci, 
6. L., Ula§, B. & Kose, O. 2007, AJ, 134, 642 

Koch, R. H., Siah, M. J. & Fanelli, M. N. 1979, 
PASP, 91, 474 

Landolt, A. U. 1964, ApJ, 140, 1494 

Landolt, A. U. 1975, PASP, 87, 409 

Lanz, T. & Hubeny, I., 2003, ApJS, 146, 417 

Massey, P., Johnson, K. E. & Dcgioia-Eastwood, 
K. 1995, ApJ, 454, 151 

Mayer, P. 1980, BAICz, 31, 292 

Mayer, P., Wolf, M., Muminovic, M. & Stupar, M. 
1986, IBVS, 2942, 1 

Mayer, P., Hadrava, P., Harmanec, P. & Chochol, 
D. 1991, BAICz, 42, 230 

Mayer, P., Niarchos, P. G., Lorenz, R., Wolf, M. 
& Christie, G. 1998, A&AS, 130, 311 

Mayer, P., Lorenz, R. & Drechsel, H. 2002, A&A, 
388, 268 

Morgenroth, O. 1935, AN, 255, 425 
Nagai, K. 2004, Var.Star Bull.Japan, 42, 1 
Paschke, A. 2010, OEJV, 130, 1 



Pearce, J. A. 1952, PASP, 64, 219 

Peter, H. 1992, BBSAG Bull., 102 

Peter, H. 1994, BBSAG Bull., 107 

Petrov, A.A. 1946, Pcrem. Zvezdy, 6, 72 

Popper, D. M. 1978, ApJ, 220, Lll 

Popper, D. M. 1980, ARA&A, 18, 115 

Popper, D. M., & Hill G. 1991, AJ, 101, 600 

Prsa, A., & Zwitter, T. 2005, ApJ, 628, 426 

Qian, S.-B., Yuan, J.-Z., Liu, L., He, J. -J. & 
Fernandez Lajiis, E., & Kreiner, J. Z. 2007, 
MNRAS, 380, 1599 

Rucinski, S. M. 1969, Acta Astron. 19, 245 

Uyamker, B., Fiirst, E., Reich, W., Aschenbach, 
B., & Wielebinski, R. 2001, A&A, 371, 675 

van Hamme, W. 1993, AJ, 106, 2096 

von Zcipel, H. 1924, MNRAS, 84, 665 

Wilson, R.E., & Devinney, E.J. 1971, ApJ, 166, 
605 

Wilson, R.E. 1979, ApJ, 234, 1054 

Wilson, R.E. 1990, ApJ, 356, 613 

Yakut, K., & Eggleton, P. P. 2005, ApJ, 629, 1055 

Yakut, K., ct al., 2013, in preparation 

Zasche, P., Uhlar, R., Kucakova, H., & Svoboda, 
P. 2011, IBVS, 6007, 1 

This 2-column preprint was prepared with the AAS IATjrjX 
macros v5.2. 



7 



Epoch (x1000) 



log T 




(b) 



4 6 8 

Epoch (x1000) 



t • • * *y • * 



4 6 8 10 

Epoch (x1000) 



Fig. 2. — (a) Residuals for the times of minimum 
light of V382 Cyg. The solid line is obtained with 
the quadratic terms and the dotted line is obtained 
with parabolic+sinus terms in ephemeris Eq. 1. 
(b) The difference between the observations and 
the quadratic ephemeris and (c) final residuals. 




Fig. 3. — Non-conservative evolution of V382 Cyg 
in the H —R diagram (top panel) and log Rvs. M 
plane (button panel). More massive and less mas- 
sive star are shown in red and green; their respec- 
tive Roche-lobe radii are dark blue and light blue. 
Initial parameters are 28.0 M©, 23.5 M and 1.72 
days with assumed Solar composition. The origi- 
nal primary is the secondary currently. 
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